The effects of spin-orbit coupling on optical properties of monolayer \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\text {MoS}_{2}$$\end{document}MoS2 due to mechanical strains

We have studied the optical conductivity of a quasi two-dimensional \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\text {MoS}_{2}$$\end{document}MoS2 in the presence of external magnetic field and spin-orbit coupling. Specially, we address the frequency dependence of optical conductivity due to spin-orbit interaction. Using linear response theory the behavior of optical conductivity has been obtained within Green’s function method. We have also considered the effects of uniaxial and biaxial in-plane strain on the optical absorption of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\text {MoS}_{2}$$\end{document}MoS2 layer. In the absence of external magnetic field with negative uniaxial strain parameter, optical conductivity includes Drude weight at zero frequency limit while Drude weight vanishes for \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\text {MoS}_{2}$$\end{document}MoS2 layer under positive uniaxial strain. Our results show that the increase of uniaxial positive strain parameter causes to move the position peak to the higher frequencies. In contrast to uniaxial strain case, the Drude weight in optical conductivity appears at positive biaxial strain value 0.15. Also we have studied the effects of magnetic field, electron doping, hole doping in the presence of spin-orbit coupling on frequency dependence of optical conductivity of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\text {MoS}_{2}$$\end{document}MoS2 in details. The magnetic field dependence of optical absorption shows a monotonic decreasing behavior for each value of temperature in the absence of strain parameter.

spin-orbit coupling of mechanical deformation has been explored in conventional GaAs based semiconductors and quantum wells where a linear strain dependence is found in the presence of this coupling 49,50 .
Strain effects have been studied to other 2D crystals and recently the possibility to tune the band gap with strain has been experimentally proven for MoS 2 51,52 . The piezoelectricity phenomenon have been investigated for monolayer MoS 2 53 so that it seems to require a deeper understanding of the effect of external nonuniform strain on the physical properties of semiconductor TMDs.
The purpose of this paper is to provide a Slater Koster tight binding model 26,48 including intrinsic spin-orbit interaction for studying the optical properties of MoS 2 monolayer in the presence of magnetic field perpendicular to the plane. Firstly, we introduce the Hamiltonian of the monolayer MoS 2 by tight-binding and then the Green's function has been obtained using the calculated energy levels. We will examine the results of optical absorption under electron/hole doping, external magnetic field, strain, and temperature. The effects of homogenous strains on the optical properties of MoS 2 plane are investigated using linear response theory. In other words, the effects of these strains on frequency dependence of optical absorption of MoS 2 layer are investigated in details. Using the suitable hopping integral and strain parameters, the band dispersion of electrons on the structure has been calculated. For calculating the optical absorption we have exploited the linear response theory in the context of Kubo formula. Green's function approach has been exploited to obtain the optical absorption of MoS 2 monolayer. The effects of external magnetic field, biaxial strains and spin-orbit coupling strength on the frequency behavior of optical absorption of MoS 2 monolayer have been studied. Also we discuss and analyze to show how spin orbit coupling, longitudinal magnetic field and biaxial strain values affect the photon frequency behavior of optical absorption of MoS 2 plane. Finally the dependence of optical absorption on biaxial strain values and spin-orbit coupling is addressed.
The parts remained of this article is planned as follows. In "Model Hamiltonian and formalism" section, we determine the tight-binding model and the Green function for the MoS 2 monolayer in the presence of the Zeeman effect, spin-orbital coupling and biaxial strain. We express our calculation of the optical absorption in Section "Optical absorption of single layer MoS 2 ". Numerical results of the optical absorption under the influence of various factors are present in Discussion and conclusion. In Section "Discussion and conclusion", includes a summary of the results and discussions.

Model Hamiltonian and formalism
MoS 2 monolayer consists of one layer of Mo atoms surrounded by two layers of S atoms in such a way that each Mo atom is coordinated by six S atoms in a trigonal prismatic geometry and each S atom is coordinated by three Mo atoms. The side and top views of the lattice structure where Mo atom is surrounded by six S atoms unit cell of MoS 2 has been indicated in Fig. 1a also we have sketched a top view of the crystal structure of MoS 2 in Fig. 1b.
In order to investigate the effect of homogenous strain and spin-orbit coupling on the electronic band structure and optical properties of monolayer MoS 2 , we apply a Slater Koster tight binding model 54 for describing the electron dynamics in the structure. This tight binding model Hamiltonian contains the relevant orbital character in the valence and conduction band. Also this model apllied for the single layer MoS 2 uses an orbital basis which includes d 3z 2 −r 2 , d xy and d x 2 −y 2 orbitals of the Mo, and p x , p y and p z orbitals of the atom S. In order to clarify the symmetric and antisymmetric hybridization of p-orbitals of S atoms localized on up and down layers, we use the definitions . The relevant physics of monolayer MoS 2 around the gap is covered by a smaller subspace, which can be obtained by performing an appropriate unitary transformation that transform the P orbitals of the top and bottom S layers into their symmetric and antisymmetric combinations with respect to the z axis. For the single-layer case, the resulting 11-band model can be decoupled in six bands with even symmetry under inversion transformation z to −z , and five bands with odd symmetry 26,55,56 . Therefore, the conduction band minimum is mainly formed from orbitals d 3z 2 −r 2 of Mo atoms and the valence band maximum is constructed from orbitals d x 2 −y 2 and d xy of  www.nature.com/scientificreports/ Mo atoms with mixing of orbitals p x , p y from S atoms 57,58 in both cases. Therefore we can consider the following basis orbitals for making Hilbert space as The tight-binding model for MoS 2 monolayer in second quantizaton representation is given by where a † i,α creates an electron in the atomic orbital α = p Sym x , p Sym y , p Anti z of S atom located in the ith unit cell. Also b † i,β implies the creation operator of electron in the atomic orbital β = d 3z 2 −r 2 , d x 2 −y 2 , d xy of Mo atom located in the ith unit cell. t Mo−Mo ij,ββ ′ introduces the hopping amplitude of electron between two Mo atoms that are located in next nearest neighbor lattice sites with atomic orbitals β, β ′ . t S−S ij,αα ′ introduces the hopping amplitude of electron between two S atoms that are located in next nearest neighbor lattice sites with atomic orbitals α, α ′ . t Mo−S ij,αβ denotes hopping amplitude of electron between Mo and S atoms that are located in nearest neighbor lattice sites with atomic orbitals α, β . Also we have introduced i, j and i, j for nearest neighbor and next nearest neighbor lattice sites respectively. ǫ S i,α refers to on-site energy of electron with atomic orbital α of S atom located in i th lattice site. In other hand, on-site energy of electron with atomic orbital β of Mo atom located in i th lattice site is presented by ǫ Mo i,β . The parameters ǫ Mo i,β in Eq. 2) with different atomic orbital basis, i.e. β = d 3z 2 −r 2 , d x 2 −y 2 , d xy are the matrix elements of a diagonal 3 × 3 matrix, ǫ Mo i . This diagonal matrix is given by Also the parameters ǫ S i,α in Eq. 2) make the following matrix presentation for ǫ S i where t ⊥ xx = V ppπ , t ⊥ yy = V ppπ , t ⊥ zz = V ppσ are the perpendicular hopping amplitudes of electron between S atoms on two different layers. A note is in order here. When we express p orbitals in terms of symmetric and antisymmetric forms, i.e. p S x , p S y , p A z , , Hamiltonian matrix elements of S orbitals will be different from the 11-band model. Under these transformations that we did, the numerical values of V ppπ and V pdπ , which are calculated within the density functional theory method, do not change 59,60 . Now we define the suitable expression for each hopping parameter in Eq. 2). All parameters in Eqs. 3,4) is given in Tables 1, 2. According to Fig. 1b, we have three nearest neighbor unit cells with connecting vectors δ l=1,2,3 ≡ R j l − R i so that R i is the position vector of ith unit cell and R j l denotes the position vector of lthe nearest neighbor unit cell. In a similar interpretation, we have three next nearest neighbor unit cells with connecting vectors b l=1,2,3 ≡ R j l − R i so that R i is the position vector of ith unit cell and R j l denotes the position vector of lthe next nearest neighbor unit cell. Depending the position vector of j l with l = 1, 2, 3 , we have three different expressions for t Mo−S ij l ,αβ in Eq. 2) as i, j l are nearest neighbor lattice sites. The matrix presentation of t Mo−S δ l with elements of t Mo−S ij l ,αβ ≡ t Mo−S δ l ,αβ for different atomic orbital basis α, β is given by (2)  Table. 2. In other hand, the position index i, j l in the matrix elements t Mo−Mo ij l ,ββ ′ (according to Eq. 2)) are related to the next nearest neighbor lattice sites. The matrix presentation of t Mo−Mo The hopping amplitudes and on-site energies have been calculated in the context of Slater-Koster scheme 26,48,52 and the values of the parameters are given in Tables 1, 2. Local spin-orbit interaction can be as well included in a suitable way 48 . The large spin-orbit in MoS 2 can be approximately understood by intra atomic contribution H SO ∝ S · L . The final results for nonzero matrix elements of H SO in Hilbert space introduced in Eq. 1) get the following relations where s = ± indicates the spin angular momentom quantum number of electrons, Mo = 86 meV and S = 52 meV imply the spin-orbit coupling strength of electron for Mo and S atoms respectively. In Eq. 8), |i, β, s� with β = d x 2 −y 2 , d xy describes the electron quantum state located on Mo atom in lattice site i with spin s in the atomic orbital labeled by β . Moreover |i, α, s� with α = p Sym y , p Sym x describes the electron quantum state located on S atom in lattice site i with spin s in the atomic orbital labeled by α . It is worthwile to some comments about only coupling between d x 2 −y 2 and d xy in Mo, and p x and p y orbitals are not considered here. Also the d 3z 2 −r 2 and p z orbitals are not involved in matrix elements of H SO . This fact can be demonstrated based on the expansion of |d x 2 −y 2 � , |d xy � and |d 3z 2 −r 2 � in terms of eigenvectors of z-component of angular momentum L z , i.e. |l, m l � in which l denotes the quantum number of total angular momentum and m l = −l, −l + 1, ..., +l refers to quantum number of operator L z . We have the following expansions for Mo atoms www.nature.com/scientificreports/ Since the spin-orbit coupling model is given by the operator form H SO ∝ S · L ∝ S z L z and also the matrix elements of operator L z are obtained as �l, m l |L z |l ′ , m l ′ � = m l δ l,l ′ δ m l m l ′ , we can expect that only the non zero matrix elements of L z are �d x 2 −y 2 |L z |d xy � and �d xy |L z |d x 2 −y 2 � . In fact |d 3z 2 −r 2 � includes only |l = 2, m l = 0� in its series expansion in terms of eigenvectors of L z , spin-orbit coupling model, i.e. S z L z , has no coupling between d 3z 2 −r 2 and d xy . In a similar reason there is no coupling between d 3z 2 −r 2 and d x 2 −y 2 . In other words �d x 2 −y 2 |L z |d 3z 2 −r 2 � and �d 3z 2 −r 2 |L z |d xy � get zero value. Also we have the following expressions for S atoms Because |p z � includes only |l = 1, m l = 0� in its series expansion in terms of eigenvectors of L z , spin-orbit coupling model, i.e. S z L z , has no coupling between p z and p x . In a similar reason there is no coupling between p z and p y . In other words �p z |L z |p x � and �p z |L z |p y � get zero value.
In the presence of applied magnetic field perpendicular to the plane of MOS 2 layer with strength B, the Zeeman splitting takes place. Based on matrix form of tight binding model Hamiltonian in Eqs. (3)(4)(5)(6)(7) and nonzero matrix elements of spin-orbit coupling model Hamiltonian, the model Hamiltonian for electrons with spin quantum number s = ± can be expressed in a more compact form once written in the k space: so that g ≈ 2.0 is the gyromagnetic constant and µ B describes the Bohr magneton constant. Also in above equation, 1 implies 3 × 3 unit matrix. Moreover the nearest ( δ j ) and the next nearest ( b j ) neighbor vectors have been shown in Fig. 1a and can be written as that b 0 = 2.41A 0 denotes the interatomic distance between Mo and S atoms. b = 3.16A 0 is the binding distance for Mo-Mo and θ B = acos( 4 7 ). The Slater-Koster tight-binding method has been applied when lattice deformations, like strain, are considered. Hence the effect of strain is driven by the dependence of tight-binding amplitudes on the interatomic distance. In the present study, the varying the interatomic bond length, as a result of applied strain, leads to the strain effects. Up to the linear order, the hopping matrixes in Eqs. (5,6,7) have been modified in the presence of strain effects as 61,62 where δ i = 7 12 b with i = 1, 2, 3 is the distance in the equilibrium positions between two nearest neighbor Mo and S atoms (Mo-S bonding length), δ i with i = 1, 2, 3 describes the distance in the presence of strain. b i = b with i = 1, 2, 3 denotes the distance in the equilibrium positions between two next nearest neighbor atoms (inplane Mo-Mo and S-S bonding lengths). Moreover b i is the distance between next nearest neighbor atoms in the presence strain. According to the definition of the electron-phonon coupling constant, the parameters β 1 , β 2 , β 3 are obtained for pp, pd, dd hybridization, respectively. The values of these parameters are β 1 = 3, β 2 = 4, β 3 = 5 . In order to study the effects of strain on the electronic properties of MoS 2 , it is necessary to express the deformed (10) www.nature.com/scientificreports/ bond lengths in terms of strain parameter and equilibrium bond lengths. We assume a general form of homogenous deformation with a large wavelength inhomogeneous strain. This strain can be considered as a powerful tool to locally manipulate the electronic structure of two-dimensional electronic systems even in atomic scale.
To consider such deformation, we use the following relations for the bond lengths, i.e. b i with i = 1, 2, 3 and δ i with i = 1, 2, 3 For homogenous biaxial strain we consider ε xx = ε yy = ε . Under applying the uniaxial strain to the MoS 2 plane, we have ε xx = −0.125ε and ε yy = ε so that ε describes the strain parameter 54 . Using the relations in Eq. 14), the deformed bond lengths, i.e. δ i and b i , are expressed in terms of strain parameter ε and undeformed bond length b. Here ω n = (2n + 1)πk B T denotes the fermionic Matsubara frequency in which T is equilibrium temperature. Total electronic density of states of MoS 2 due biaxial strains and under applying external magnetic field can be obtained by electronic band structure as Summation over wave vectors have been performed into first Brillouin zone of honeycomb lattice. The density of states includes prominent asymmetric peaks due to the band edge of parabolic subbands. The peaks positions arise from the band edge state energies and the density of states heights are proportional to inverse square root of the subband curvature and band degeneracy. For determining the chemical potential, µ , we use the relation between concentration of electrons ( n e ) and chemical potential. This relation is given by Based on the values of electronic concentration n e , the chemical potential, µ , can be obtained by means Eq. (19).

Optical absorption of single layer MoS 2
The optical conductivity is obtained as the response of the electrical current ( J e ) to an external electrical field. Imposing the continuity equation for the charge density ρ , i.e. ∂ ∂t ρ + ∇ · J e = 0 , the explicit form of the electrical current operator can be calculated. This calculation has been done for a bilinear Hamiltonian describing Fermionic system 63 and we can exploit this result for model Hamiltonian in Eq. (15). The operator form of electrical current operator J e for itinerant electrons of MoS 2 layer is given by (14)  www.nature.com/scientificreports/ so that v s,η (k) = ∇ k E s η (k) denotes the group velocity of electrons with spin s at band structure η . The linear response theory is implemented to obtain the optical conductivity under the assumption of a low dynamical electric field (as a perturbing field). The Kubo formula gives the transport coefficient σ (ω) in terms of a correlation function of electrical current operators where it is assumed that electrical current flows along zigzag direction, i.e. x direction in Fig. 1,. By implementing Wick's theorem, we can calculate the correlation function between current operators in Eq. (21) as Applying the Wick's theorem leads to the following expression for energy current correlation function as

Discussion and conclusion
In this section, we have presented the numerical results of optical absorption rate of electromagnetic waves in MoS 2 monolayer in the presence of magnetic field, spin-orbit coupling and homogenous strain effects. We investigate the frequency dependence of optical conductivity Im(ǫ(ω)) of MoS 2 layer due to variation of physical parameters. Using dimensionless homogenous strain parameter ε , the amounts of redefined hopping amplitudes of MoS 2 layer, i.e. t Mo−S   www.nature.com/scientificreports/ Using the electronic band structure of MoS 2 monolayer, we can obtain the electronic Green's function in Eq. (17). It should be noted that the variation of chemical potential µ leads to various electronic concentration according to Eq. (19). Afterwards density of states, dynamical electrical conductivity σ (ω) are found by substitution of Green's function into Eq. (26) so that optical absorption rate can be found by Eq. (27). In our numerical results we assume the polarization of electric field of electromagnetic wave is along the zigzag direction according to Fig. 1. Both inter and intra band transitions contribute to the results of optical properties of MoS 2 monolayer. The distance between two next nearest neighbor atoms in the atomic structure of the MoS 2 layer is assumed to be b = 1.
The frequency dependence of optical absorption of single layer undoped MoS 2 in the absence of strain for different values of magnetic field has been shown in Fig. 2. The temperature has been fixed at amount k B T = 0.1 eV. Based on this figure, it is clearly observed that the Drude weight value increases with magnetic field so that the Drude weight for gµ B B = 0.3 eV gets around normalized value 15. The origin of this Drude weight or the zero frequency limit of optical absorption comes from intraband transition of electrons due to classical behavior of them in this limit. In other hand the peak in optical absorption for all amounts of magnetic field appears at finite frequency ω ≈ 2.5 eV. The height of this peak decreases with magnetic field. Such peak arises from interband transition of electrons. An additional peak appears in optical conductivity for magnetic fields gµ B B = 0.2, 0.3 eV. This additional peak is located at higher frequencies rather than first peak position. The frequency position of second peak in optical conductivity for gµ B B = 0.2 eV ( gµ B B = 0.3 eV) is found around ω ≈ 2.85 eV ( ω ≈ 3.1 eV). Also the optical absorption decreases with frequency in the region ω > 3.25 eV. Figure 3 shows the dependence of optical absorption Im(ǫ) on magnetic field gµ B B for different values of temperature k B T at fixed frequency ω = 2.5 eV in the absence of any type of strain, i.e. ε = 0.0 . A monotonic decreasing behavior for magnetic field dependence of optical absorption for each value of temperature is clearly observed. It can be understood from this point that increasing magnetic field leads to enhance the band gap in density of states and consequently optical absorption of single layer MoS 2 decreases with magnetic field. At fixed magnetic field in the region gµ B B > 0.8 eV, the increase of temperature leads to enhance the transition rate of electrons so that optical absorption rises in this magnetic field region. Also Fig. 3 implies at fixed magnetic field in the region gµ B B < 0.8 eV, the increase of temperature reduces the optical absorption of MoS 2 layer due to scattering rate of electrons. Moreover the optical absorption curves for k B T = 0.05 eV, 0.07 eV, 0.1 eV fall on each other at magnetic fields above 2.0 eV.
The behavior of chemical potential dependence of optical conductivity of doped MoS 2 layer for different photon frequencies, namely ω = 0.5, 1.0, 1.5, 2.0 , has been shown in Fig. 4. The temperature has been fixed at k B T = 0.1 eV and any type strain parameter value is assumed to be zero. This figure implies that there is a peak in chemical potential dependence of Imε(ω) for frequency value ω = 0.5 . At chemical potential values µ > 0.8 eV for ω = 0.5 , the increase of chemical potential leads to enhance the electronic concentration so that the scattering rate between electrons rises and consequently Im(ǫ) reduces with µ in this region. In other hand at chemical potentials below 0.8 eV, the increase of µ causes to the transition rate of electrons from ground state to excited ones. This fact denotes the increase of optical absorption with chemical potential in the region µ < 0.8 at ω = 0.5 . However there is no peak in optical absorption for the other frequency values, i.e. ω = 1.0, 1.5, 2.0.
We have also studied the effects of temperature on chemical potential dependence of Im(ǫ(ω)) . In Fig. 5, we have presented the behavior of optical absorption of doped MoS 2 layer in terms of µ in the absence of any www.nature.com/scientificreports/ strain parameter for different values of temperature. The photon frequency has been fixed at ω = 2.5 eV. This plot indicates a monotonic decreasing behavior for each value of temperature. It can be understood from the fact that the increase of chemical potential leads to enhance the electronic concentration. Thus the scattering rate between electrons rises which reduces the optical conductivity. Another feature is pronounced in this figure. For fixed chemical potential below 0.8 eV, Im(ǫ(ω)) decreases with temperature according to Fig. 5. Higher temperatures cause more scattering of electrons which reduces the optical conductivity. However for chemical potentials above 0.8 eV, the increase of temperature leads to enhance the transition rate of electrons between quantum energy levels. This fact implies optical conductivity increases with temperature in the chemical potential region above 0.8 eV.
We have studied the effect of in-plane uniaxial strain along armchair direction with strain parameter ε on frequency dependence of Im(ǫ) for undoped single layer MoS 2 structure in Fig. 6. In this figure, the effects of  The applied longitudinal magnetic field is assumed to be zero and temperature has been fixed at k B T = 0.1 eV. In the absence of strain parameter, a peak appears at finite frequency ω = 2.4 eV due to interband transition effects of electrons. For finite non zero strain parameter, two peaks in optical absorption appear at finite frequencies. According to Fig. 6, the intensity of optical absorption at peak frequency position decreases with uniaxial strain. Also the distance between two peaks in Im(ǫ) at finite uniaxial strain parameter increases with ε as shown in Fig. 6. Moreover the figure indicates that the optical absorption vanishes in frequency region ω < 1.5 for all values of ε . Also there is no Drude weight at zero frequency limit of optical absorption for all uniaxial strains due to band gap in excitation spectrum of the electronic band structure. In fact, the intraband transition in the presence of uniaxial strain has no considerable contribution to the electronic   Fig. 7. The temperature and magnetic field have been assumed to be k B T = 0.1 eV and gµ B B = 0.0 eV. There are two peaks in optical absorption at finite frequency for each non zero strain parameter so that the distance between peaks increases with absolute value of ε . Only for ε = −0.15 , optical absorption gets the non zero value at zero limit frequency. In other words, there is non zero Drude weight in optical absorption for ε = −0.15 . Such fact demonstrates that intraband electronic transition contributes to the optical absorption and consequently the system behaves as a metal. MoS 2 layer under negative biaxial strain behaves as a transparent media at frequencies ω > 4.0 eV so that optical absorption gets the zero value in this frequency region. Also Im(ǫ) vanishes for ε = 0.0, −0.05, −0.1 at frequencies ω < 0.7 eV. www.nature.com/scientificreports/ The effects of in-plane biaxial strain ε xx = ε yy = ε with positive sign on behaviors of optical absorption have been studied in Fig. 8. We have plotted the frequency dependence of Im(ǫ(ω)) for different values of positive in-plane biaxial strain ε in the absence of magnetic field and temperature k B T = 0.1 eV in Fig. 8. The intraband band transition contributes to the optical absorption at zero frequency limit for ε = 0.15 so that Drude weight takes the remarkable value for this strain amount. In addition to the Drude weight in optical absorption, the finite frequency peaks in optical absorption curves are clearly observed. Such finite frequency peaks in optical absorption arises from electronic interband transitions. The frequency position of peak in optical conductivity moves to lower amounts with ε . This can be understood from this fact that the increase of ε causes the decrease of band gap in density of states. Thus the peak in Im(ǫ(ω)) appears at lower frequency with increase of strain parameter ε . Another novel feature in Im(ǫ(ω)) is pronounced in Fig. 8. The optical absorption gets the zero value at frequencies ω < 1.0 eV for ε = 0.0, 0.05 . The metallic property of MoS 2 layer in the presence of strain value ε = 0.15 is considerable since the Drude weight in optical conductivity for strain parameter ε = 0.15 gets higher value in comparison with the other values of biaxial strain parameters.
The behaviors of frequency dependence of optical absorption of undoped MoS 2 layer for different negative biaxial strains have been presented in Fig. 9. The temperature and magnetic field have been assumed to be k B T = 0.1 eV and gµ B B = 0.0 eV. There is a finite frequency peak in optical absorption for each value of strain parameter so that the peak position tends to higher frequency with absolute value of ε . According to Fig. 9, Drude weight of MoS 2 layer gets a very low value for all amounts in-plane biaxial strain. It can be justified from this point that intraband electronic transition gives no considerable contribution to the optical absorption and consequently the system behaves as a insulator. MoS 2 layer under negative biaxial strain behaves as a transparent media at frequencies ω > 4.0 eV and ω < 1.5 since Im(ǫ) gets the zero value in these frequency regions.
We have studied the hole doping effects on the behavior optical absorption of MoS 2 . For this purpose, the dependence of Im(ǫ(ω)) on the negative chemical potential µ for different frequencies ω at fixed temperature k B T = 0.1 eV in the absence of magnetic field has been plotted in Fig. 10. A considerable peak in optical conductivity MoS 2 layer is clearly observed for ω = 0.5 eV. The chemical potential dependence for optical conductivity at ω = 2.0 eV indicates that Im(ǫ(ω)) shows a monotonic increasing behavior with decrease of absolute value in the region µ < −1.0 . Upon more increase of µ above −1.0 , optical conductivity for ω = 2.0 eV reaches a constant value according to Fig. 10.
In Fig. 11, we have presented the behavior of optical absorption of doped MoS 2 layer in terms of negative chemical µ in the absence of any strain parameter for different values of temperature. The photon has been fixed at ω = 2.5 eV. This plot indicates a monotonic increasing behavior for each value of temperature. It can be understood from the fact that the decrease of absolute value of chemical potential leads to reduce the hole concentration. Thus the scattering rate between holes decreases which increases the optical conductivity. Another feature is pronounced in this figure. For fixed chemical potential above −1.6 eV, Im(ǫ(ω)) decreases with temperature according to Fig. 5. Higher temperature causes more scattering of holes which reduces the optical conductivity. However for chemical potentials below -1.6 eV, the increase of temperature leads to enhance the transition rate of electrons from valence band to conduction one. This fact implies optical conductivity increases with temperature in the chemical potential region below −1.6 eV.
Finally we have plotted the electronic band structure of MoS 2 monolayer in the absence strain parameters for various external magnetic fields in Fig. 12. This figure implies the increase of applied magnetic field leads  www.nature.com/scientificreports/ to decrease distance between energy levels as the structure tends to metallic property with magnetic field. Such metallic property with magnetic field has been approved in Fig. 2 where the Drude weight in optical conductivity increases with magnetic field.

Data availability
The datasets used and analysed during the current study available from the corresponding author on reasonable request.